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ABSTRACT 

We study the radial distribution of supernova la (SNe la) in morphologically 
selected early-type host galaxies from the Sloan Digital Sky Survey (SDSS) and discuss 
its implications for the progenitor systems of SNe la. While new observations of early- 
type galaxies suggest that they contain small fractions of young stellar populations, 
they are also the most likely hosts for long time delay SNe la. We find that there 
is no statistically significant difference between the radial distribution of SNe la and 
the light profile of their early-type host galaxies, which are dominated by old, metal- 
rich stellar populations. This confirms the commonly accepted idea that some SN la 
progenitors have time delays of the order of several Gyr. 
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1 INTRODUCTION 

Type la supernovae (SNe la) have been exte nsively used 
as di s tance indicators for co smology (see e.g. IRiess et al.l 
Il998l : iPerlmutter et al.lll999l ) using the empirical Phillips 
relation i|Phillipslll993i ) and related methods. They are also 
important for understanding galaxy evolution, affecting the 
chemistry and energy budget of their hosts, and the star 
formation and metallicity evolution of the Universe. Uncer- 
tainties in both the progenit or scenarios and the explosion 
mechanism (for a review see iHillebrandt fc Niemeverll2000l ) 
prevent us from having a convincing physical picture that 
justifies their use as standardisable candles at high redshift 
and is able to reproduce their observed diversity in the local 
Universe. 



1.1 The SN la time delay distribution 

Linking SN la progenitor scenarios with cosmology, galaxy 
evolution and the cosmic star formation history requires in- 
formation on the time between their formation and explo- 
sion, or time delay distribution. Thermonuclear explosions 
seem to originate in white dwarf stars (WDs), with main 
sequence (MS) lifetimes ranging from ~ 30 Myr to several 
billion years, thus the minimum time delay must be of the 
order of ~ 30 Myr. To narrow the constraints details of pos- 
sible progenitor systems need to be used. 
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Most proposed progenitor scenarios involve mass trans- 
fer on to a CO WD that reaches the Chandrasekhar limit in 
a binary system. In these scenarios the WD grows through 
either the expansion and Roche lobe overflow of an evolved 
companion {single degenerate [SD] scenarios), or the slow 
release of gravitational waves, orbital shrinking, Roche lobe 
overflow and merging of a compact double WD system 
(double-degenerate [DDJ scenarios). 

For the SD scenario it has been suggested that 
only systems that can transfer enough mass under a 
critical rate ca n bring the WD star to the Chan- 
drasekhar limit (|Nomoto fc Kondol Il99li ). In binary pop- 
ulation synthesis (BPS) simulations most of the WDs 
that reach the Chandrasekhar mass accrete matter from 
a slightly evolved MS star, the so called CO WD + 
MS - SD scenario (lHachisu et al.lll999l : lLanger et a"ill2000l : 
lHan fc Podsiadlowskill2004l ). In this channel the mass of the 
companion determines both when accretion starts and the 
rate of accretion. As a consequence, the time-delay distribu- 
tion of this channel seem to be relatively narrow, peaking at 
~ 670 Myr and rapidly becoming negligible after ~ 1.5 Gyr. 

Several Gyr time-delay progenitors can also be pro- 
duced in the SD scenario when a relatively low mass red- 
giant (RG) star transfers matter on to a CO WD star, i n 
the CO WD + RG - SD scenario |Hachisu fc Nomotdll996l ). 
This s cenario appears to be supported by observational evi- 
dence (IPatat et al.ll2007|). although its relative contribution 
is unclear ( Han fc Podsiadlowskill2004 ). 

The characteristic time delay in the DD scenario 
l|lben fc Tutukovl 1 1984 IWebbinkl [l98i ) is the coalescence 
time-scale of the binary system, which depends roughly on 
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the fourt h power of the separa t ion o f the double-degenerate 
system (Sha piro fc Teukolskvl Il983h . The time-delay dis- 
tribution can be described by a low time-delay cutoff (~ 
30 — 100 Myr), an approximately power-law decline up to 
the age of the Universe and some secondary peaks depending 
on the details of the BPS simulation. 

However, the expected accretion rates in the DD sce- 
nario are thought to be too big to lead to successful ther- 
monuclear explosions, with most simulations leading to 
accretion-induced collapse (AIC) and the formatio n of com- 
pact objects instead (see ISaio fc Nomotd Il998l . and ref- 
erences therein). Rotation or more sophisticated neutrino 
physics may give this channel a more significa nt role in the 
future jYoon, Podsiadlowski. fc Rosswod 120071 ) . 

Observationally, a comparison between the cosmic 
supernova rate and star formation history has led 
to indirec t constraints on t he ti m e delay distr i bution 
fsee e.g. iGal-Yam fc Maozl 12004 IStroker et all 12004 
iHopkins fc Beacoml 20061 ). but uncertainties in both func- 
tional quantities seem unlikely to give constraints with the 
required accuracy to distinguish between d ifferent channels 
|Forster et al.ll2006l ; iBlanc fc Greggidl200Sl ). 

A better probe for clues regarding the SN la time de- 
lay distribution might be to look at t heir individual host 
galaxies (see e.g 



Mannucci et alj|2005l; ISullivan et "aill2006l ; 
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IPrieto. Stanek. fc Be acom 2008) and their galaxy environ- 
ments (see e.g. Isharon et alJ libo?). With this approach un- 
certainties in the star formation histories of individual galax- 
ies can be better controlled if accurate stellar population 
reconstructions are used. 



whether the radial distribution of SNe la follows a de Vau- 
couleurs law. 



2 SAMPLE SELECTION 

We create a sample of SN la early-type host galaxies by se- 
lecting all SNe la in the CfA list of supernovaqj whose hosts 
have been covered by the Sloan Digita l Sky Survey (SDSS, 
I York et al.ll2000l : IStouehton et al.ll2002l ). The SDSS is a sur- 
vey of half the Northern Sky, providing us with photome- 
try in the five optical f ilters u,g,r,i and z dFukugita et al.l 
Il996l : iGunn et al.lll998l). We use the current data release 
DR6 (jAdelman- McCarthy et~al1l200St ). 

We visually inspect every SN la host galaxy observed 
by SDSS and select only those with early-type morphology. 
Visual classification avoids the introduction of biases asso- 
ciated with proxies for morpholo gy such as colour, c oncen- 
tration or structural parameters l)Lintott et al.ll2008l ). From 
the remaining galaxies we select those with half-light radii 
bigger than 2.4", or six pixels in the SDSS images. The final 
sample contained galaxies with a mean and median half- 
light radius of 17.6" and 12.9", or about 44 and 32 pixels, 
respectively. Their redshifts were in the range 0.001 to 0.3, 
with a mean of 0.06 and median of 0.05. 

About three quarters of the SNe in the final sample were 
dicovered by the KAIT/LOSS SN survey, the SDSS SN sur- 
vey and amateur astronomers in almost equal proportions. 
The rest were discovered by various surveys including the 
Nearby SN Factory, among others. 



1.2 Early— type galaxies 

The dominant stellar populations of early-type galaxies are 
old and must have f ormed in the high redshift universe (see 
iThomas et al.l [20051 . and references therein). 

However, the GALEX UV space telescope has led to the 
discovery that a significant fraction of early-type galaxies in 
the low redshift unive rse formed a few percent of their stellar 
mass in the last Gyr llYi et al . 2005; Sc hawinski et al.1l2006l : 
iKavirai et al.ll2007l ; ISchawinski et alj|2007l ). These episodes 
of residual star formation are not apparent in the optical, but 
become prominent in the UV. This opens up the possibility 
that the SNe la observed in apparently passive early-type 
galaxies may be due to these small young populations and 
so might have short time delays. 

The bulk stellar population of massive e arly- type galax- 
ies fo llow a r 1//4 or de Vaucouleurs law l)de Vaucouleursl 
1 19481 ) and the dynamical relaxation time-scale of these 
systems is much higher than the age of the Universe 



l)Binnev fc Tremaind 1983). High resolution H/3 maps have 
been used to trace young stellar populations in early- 
type galaxies and they tend not to follo w the typ- 
ical de Vaucouleurs profile of their hosts (ISa rzi et al.l 
120061 : iKuntschner" et alJ [2006). Evidence for central molec- 
ular disks in early- type galaxies supp orts this idea 
l|Lucero fc Youndl2007l : ICrocker et al.ll2008l ). 

If the radial distribution of SNe la follows the light of 
the bulk stellar population in early-type galaxies, then their 
progenitors must have an age of at least the minimum be- 
tween the relaxation time-scale and the age of the bulk of 
the stellar population. We attempt to quantitatively test 



3 ANALYSIS 

We use elliptical projected radial coordinates and model the 
surface bri ghtness prof iles of the early-type galaxies in our 
sample as (Ser sidll968l ): 
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where n is known as the Sersic index (normally assumed to 
be four in early-type galaxies) and by definition half of the 
light comes from r < r c , with I c = I(r a ). These requirements 
define b n , approximately 2n — 1/3. 

l/r, 

Now, using the change of variable u = 6 n ( ■&■ 1 one 



can show that the probability density per unit variable u of 
a photon coming at a radius r = r c (u/b n ) n is: 



fu(u) = 



f °° t^-ie-tdt T(2n) 



(2) 



where T(u) = J™ t^e^dt. This allows one to define an 
equivalent SN la rate surface magnitude brightness: 
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which, for a given n, should have the same slope of the log 
of the true surface brightness in (r/r c ) 1//n space, but differ 
by a constant. 

1 See: http: //www. cf a. harvard. edu/iau/lists/Supernovae .html 
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Thus, the cumulative probability distribution per unit 
variable u that a photon comes from inside a region of radius 
r = r c (u/b n ) n is: 



SNe la with de Voucouleurs profile 



F u (u) = P(2n, u), 



(4) 



where P(a, u) is the regularised incomplete Gamma func- 
tion defined as P(a, u) = f^t^e^dt / T(a). 

With the distribution of normalised distances of ob- 
served SNe la it is possible to test whether the SN la rate is 
proportional to the surface brightness intensity in early-type 
galaxies. 



3.1 Systematic effects 

Since core collapse SNe are very rare in early-type hosts, SN 
type misclassification is unlikely. However, early-type galax- 
ies have extended radial profiles that easily overlap with 
other galaxies in the line of sight, making host galaxy iden- 
tification problems important at large radii. Moreover, sev- 
eral effects will occur at small radii: different point -spread 
functions, deviations from the Sersic profile close to galaxy 
cores, astrometric errors betwe en SPSS an d SN discovery 
coordinates, or the Shaw effect (|Shawlll979r ). i.e. undetected 
SNe near the core of galaxies due to the saturation of old 
photographic plates. To minimise these systematic errors we 

limit our sample to SNe with intermediate radii. 

SDSS positions are accurate to 0.1" l|Pier et al.ll2003l ) 
and only 15% of the SNe in our sample were discovered be- 
fore 2000, which means that astrometric problems or the 
Shaw effect are unlikely to be significant. We define a min- 
imum SN distance of 0.2 half-light radii, since this central 
zone is well resolved in more than 80% of the images. For the 
smallest allowed half-light radii this corresponds to about 
five times the SDSS astrometric accuracy, but about two 
times the SDSS resolution. We also define an arbitrary max- 
imum distance of four half-light radii, which is the radius at 
which approximately 85% of the light of the galaxy is con- 
tained in a de Vaucouleurs profile. The resulting mean and 
median SN distances from the galaxy cores were 18.2" and 
9.6", about 45 and 24 pixels, or 13 and 7 times de resolution, 
respectively. 

Thus, we limit our sample to SNe with 0.2 < r/r c < 4. 
When more than one host candidate met these requirements 
we arbitrarily chose the one closer to the SNe. This could 
bias our results to a centrally concentrated distribution, but 
dim host galaxies that are not seen in SDSS images could 
have the opposite effect. 

In Fig. [T]we plot the equivalent of the surface brightness 
profile using type la SNe and assuming a de Vaucouleurs 
law with the restricted normalised distances defined above. 
As a first approximation it shows that SNe la follow their 
early-type host light profiles. Other authors have studied 
the radial distribution of type la SNe, but not using half- 
light r adii normalised distances (llvanov. Hamuv. fc Pinto! 



20071 ). Recently, 



200d: iBartunov, Tsvetkov. fc Pavlvuk 
Totani et al.l ( 20081 ) have also used half-light radii nor- 
malised distances, but with photometrically typed SNe la 
at high redshift and with a much lower relative resolution. 




(r/O 



Figure 1. Observed (error bars) values of /^j a , defined by cq. 
vs the expected values (line) if SNe la follow a de Vaucouleurs 
law in early— type galaxies. The binning is robust at intermediate 
radii, with 5 or more elements in more than 80% of the bins. 



3.2 Goodness of Fit test 

To quantitatively test the hypothesis that SNe la follow 
their early-type hosts light profiles in this limited sample we 
normalise the cumulative probability distribution in equa- 
tion (0}, so that 6 n 0.2 1/n < u < 6 n 4 1 /™, and contrast it with 
the observed cumulative distribution using a Kolmogorov- 
Smirnov (KS) goodness of fit test. We conclude that for this 
restricted region the SN la rate is statistically consistent 
with a de Vaucouleurs profile (see Fig. [2}. 

However, it is known that early-type galaxies light pro- 
files have a continuous distribution of Sersic indices connect- 
ing dwarf to giant ellipticals, increasing from approximately 
n = 2 to n = 6 (|Graham fc Guzma.nl [20031 ). We have in- 
cluded a variable Sersic index in our fitting procedure and 
found that the best-fitting Sersic indices are distributed in 
approximately the same range, but have a slight preference 
for lower values, with an average of 3.4 and a median of 3.3. 
This distribution might be a result of the way SNe la have 
historically been discovered, which in many cases depends 
on pre-selected galaxy samples. 

To allow for this diversity to be included in our analysis 
we use a new change of variable suggested by eq. @ , namely 
v = P(2n,u), or the fraction of the integrated light of the 
galaxy at the position of the SNe if the profile was well 
described by equation JTJ. The new probability density and 
cumulative probability density will be simply: 



and 



F v (v) 



(•») 



which are independent of n. This allows us to compare all 
galaxies irrespective of their Sersic indices. 

Hence, after finding the best-fitting values of n for each 
individual galaxy we set a cut for the SN sample based on 
the new variable v rather than u. Given that the half-light 
radii resulting from the Sersic profile fitting were on average 
smaller, we limit our sample to SNe within 0.2 < u < 0.85, 
equivalent to 0.3 < r/r c < 4 for the average Sersic index. 
We also select galaxies were the best-fitting Sersic indices 
were in the expected range 2 < n < 6. The distribution of 
the variable v contrasted to F v (v) is in Fig. [2] 
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Figure 2. Kolmogorov Smirnov tests of the cumulative distribution of SNe la (top) and SN la histograms (bottom) assuming a de 
Vaucouleurs profile (left) or a Sersic profile (right) and the predicted distributions. The data was binned so that more than 80% of the 
bins would have at least 5 elements. Note that different selection criteria led to different sample sizes. 



4 CONCLUSIONS AND DISCUSSION 

We have found that there is no statistically significant dif- 
ference between the radial dependence of the SN la rate in 
morphologically selected early-type galaxies and their host 
surface brightness profiles at intermediate radii. The proba- 
bilities of not rejecting the hypotheses that SNe la follow a 
de Vaucouleurs or Sersic profile is 24.0% and 27.0%. 

In this work we use host half-light radii (r e ) normalised 
distances, obtained after fitting de Vaucouleurs or Sersic 
profiles on r band SDSS images, and quantified any devia- 
tions from a de Vaucouleurs or Sersic profile at intermediate 
radii. Intermediate radii are defined in Section [3] 

The main implications of this work is that the time de- 
lay of type la SNe in early-type galaxies should be of the 
order of the age of their host galaxies, or several Gyr. This 
assumes that young stellar populations of early-type galax- 
ies do not follow the light distribution of their hosts, which is 
supported by recent observations of central molecular disks 
and high-resolution SAURON H/3 maps and stellar popu- 
lations reconstructions, showing that when young (<1 Gyr) 
stellar populations are present, they are mainly in the cen- 
tral regions of their hosts. Even after including small radii 



in our sample we did not see an excess of the SN la rate 
in the central parts of early-type galaxies, but given all the 
possible systematic errors of our method we cannot conclude 
if there is any real deficit. 

Our result does not exclude either the SD or DD sce- 
narios in a picture of multiple progenitor scenarios, but con- 
firms the idea that very young SN la progenitor scenarios are 
probably not dominant in early-type galaxies. They rather 
indicate that long time delays of several Gyrs for SNe la in 
early-type galaxies are dominant. 

Interestingly, early-type galaxies are known to ex- 
hibit radial gra dients in me t allicit y, but not in a ge or 
Q-enhancement (|Saglia et all |2000| ; iMehlert etHl 120031 ; 
Wu ct al. 2005). The metallicity gradients range all the way 
to the outer parts of early- type galaxies, beyond two effective 
radii i|Mendez et all [2005). Given this, our result also sug- 
gests that within the metallicity range exhibited by massive 
early-type galaxies, the SN la rate appears not to be strongly 
affected by metallicity, which hints that long time-delay SN 
la progenitor scenarios may not have a very strong metallic- 
ity dependence. Note that a significant dep endence has been 
predicted for metallicities below solar (see iKobavashi et al.l 
ll998l ; lMeng. Chen, fc Hanll200Sl ). 
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If the metallicity gradients and surface brightness pro- 
files are relics of the process of galaxy formation in early- 
type galaxies, our result would suggest that the progenitors 
of long time-delay SN la are as old as their host galaxies. 

Future studies of the radial distribution of sub-classes 
of SNe la, e.g. SN1991bg and SN1991T-like events, should 
provide further insights in to whether there are multiple p ro- 
genitor scenarios (see e.g. IScannapieco fe B ildstcn 2 0051 ). 
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